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Calculations of reaction rates for the third-order QED process of photon splitting 7 ^ 77 
in strong magnetic fields traditionally have employed either the efi'ective Lagrangian method or 
variants of Schwinger's proper-time technique. Recently, Mentzel, Berg and Wunner [Q] presented 
an alternative derivation via an S-matrix formulation in the Landau representation. Advantages of 
such a formulation include the ability to compute rates near pair resonances above pair threshold. 
This paper presents new developments of the Landau representation formalism as applied to photon 
splitting, providing significant advances beyond the work of Q by summing over the spin quantum 
numbers of the electron propagators, and analytically integrating over the component of momentum 
of the intermediate states that is parallel to field. The ensuing tractable expressions for the scattering 
amplitudes are satisfyingly compact, and of an appearance familiar to S-matrix theory applications. 
Such developments can facilitate numerical computations of splitting considerably both below and 
above pair threshold. Specializations to two regimes of interest are obtained, namely the limit 
of highly supercritical fields and the domain where photon energies are far inferior to that for 
the threshold of single-photon pair creation. In particular, for the first time the low-frequency 
amplitudes are simply expressed in terms of the Gamma function, its integral and its derivatives. 
In addition, the equivalence of the asymptotic forms in these two domains to extant results from 
effective Lagrangian/proper-time formulations is demonstrated. 
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I. INTRODUCTION 



The third-order quantum electrodynamical process of photon splitting 7 — s- 77 in a strong magnetic field, currently 
popular in several astrophysical models of different neutron star sources, was first studied over three decades ago. 
Due to analytic complexities encountered when investigating this interaction, it was not until the beginning of the 
1970s that a body of correct and uncontroversial results emerged. These early splitting calculations used either 
effective Lagrangian or variations of Schwinger's proper-time techniques [^-0, the expediency of which yielded 
compact analytic forms for the rates R when specializing to low energy ( i? oc ) or low field ( R (x B^) cases. 
After a hiatus of nearly two decades, photon splitting became of interest again in the literature following the 

publication of an S-matrix calculation in the Landau representation of its rates by Mentzel, Berg and Wunner [Q, 
specifically because of their contention that the earlier works cited above had seriously underestimated the strength 
of this process. The rates computed in Q were later retracted in p^, with a sign error in their numerical coding 
having been discovered and corrected. Mentzel et al.'s analytic derivation was the first comprehensive presentation 
of the application of a Landau representation technique specifically to magnetic photon splitting, though the QED 
formalism presented by Melrose and Parle 1 13 1^ virtually provided an equivalent enunciation of such S-matrix forms 
for splitting amplitudes. More recently, Weise, Baring & Melrose [|ll[ confirmed the analytic derivation of |]l|. The 
Landau representation calculations and most of the earlier effective Lagrangian and proper-time presentations were 
generally applicable to non-dispersive regimes below the pair creation threshold ( hu! — 2mc^ ) , where the momentum 
vectors of the initial and final photons are coUinear, and arbitrary field strengths. 

Below pair threshold, the effective Lagrangian approach of and the proper-time calculations in 1^-^ appear 
much more amenable for the purposes of numerical evaluation than the S-matrix formulation in the Landau repre- 
sentation. This arises because effective Lagrangian and proper-time (collectively referred to by the label ELP here) 
methods produce results that involve triple integrals over relatively simple (hyperbolic and exponential) functions, 
while the S-matrix amplitudes integrate over the parallel momentum pz and include a triple summation over the 
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Landau level quantum numbers of the intermediate pair states. Both techniques start from different but equivalent 
forms of the electron propagator, and hence S-matrix computations should yield identical results to proper- 

time numerics [^|-^. For the specific case of magnetic pair creation 7 e+e^ , such an equivalence of the S-matrix 
and proper-time methods has been demonstrated but only via continuous asymptotic approximations that 

smoothly average out the exact "sawtooth" resonance structure. Yet the S-matrix Landau representation approach 
explicitly retains the resonances in the scattering amplitudes above pair threshold, whereas the ELP methods elimi- 
nate such information early during developments. Photon splitting becomes effectively first-order in af at any one of 
a multitude of pair resonances, generated when the intermediate states become "on-shell." Hence it is quite possible 
that splitting can compete effectively with pair creation as a photon absorption mechanism above pair threshold. As- 
certaining whether this is true is an interesting physics question. Moreover, if splitting is approximately as probable 
as pair creation above threshold, then it manifestly changes the character of vacuum dispersion, so that quadratic 
(and by inference perhaps higher order) contributions to the vacuum polarization tensor become significant relative 
to the standard linear ones used in the derivation Q of kinematic selection rules for splitting. Hence, the generation 
of exact and compact expressions for the rates for 7 — > 77 valid both below and above pair creation threshold is 
clearly a worthwhile enterprise from a physics perspective. 

Developed expressions for the rates for photon splitting are also important for astrophysical applications of this 
process, particularly to effect efficient and accurate computations of such rates. These applications have so far focused 
on neutron star magnetospheres, primarily on models of soft gamma repeaters (SGRs) and strongly-magnetized 
pulsars, both being extremely topical in the astrophysics community at present. The potential importance of splitting 
in neutron star environments was suggested by [p| Jl^ , p^ . Possible formation of splitting cascades has been explored 
in models of SGR transient outbursts as a means of softening the spectrum efficiently with no production of pairs 
p0|-p3|. If both polarizations can split, or if polarization switching is active during SGR outbursts, then the properties 
of the splitting cross-section guarantee emergent spectra in the observed range (20-150 keV) and of the observed shape 
for all fields in excess of around 10^** Gauss pO|-p^, provided that the emission region is not concentrated near the 
polar cap. The spectral properties of SGRs in quiescent emission appear to be distinct from those during outburst. 
Pulsations and temporal increases of their periods (i.e. spin-down) have now been observed |24| - [2^ for two of the four 
confirmed SGRs (SGR 1806-20 and SGR 1900-^14), leading to inferences of fields in the vicinity of lO^^ Gauss. The 
connection between these pulsars of extremely high magnetization, so-called magnetars, and conventional radio/X- 
ray/gamma-ray pulsars is not well-understood. Baring & Harding [ p7[ postulated that radio quiescence, a property of 
the SGRs, may be common in magnetars due to the efficient action of photon splitting and other effects in suppressing 
the creation of pairs. Photon splitting also has spectral implications for such pulsars with more modest fields: |2|] 
demonstrated that the unusual absence of > 30 MeV emission in the gamma-ray pulsar PSR 1509-58 (whose spin- 
down field is ~ 3 X 10^"^ Gauss) can naturally be explained by the operation of 7 ^ 77 in the intense magnetic and 
gravitational fields near its surface. 

Several desirable goals are immediately identifiable on the basis of this historical path for the study of the physics of 
photon splitting, and the needs of the astrophysics community. It would be satisfying (i) to obtain analytic expressions 
for rates that are valid above pair creation threshold using the Landau representation methodology, (ii) to know 
whether the analytic formalism of Mentzel et al. |Q| can be developed and simplified, and (iii) to demonstrate a formal 
equivalence between this S-matrix Landau representation approach and extant results from proper-time/effective 
Lagrangian techniques. This paper addresses these issues, using the verified analytic formalism of Mentzel et al. as 
the starting point for mathematical developments. The analysis here considers all the polarization modes that are 
permitted by the CP invariance symmetry (_L— !-|||| , _L^_L_L and ||— >-L|| ), and applies for coUinear momenta of the 
incoming and outgoing photons, i.e. when the effects of vacuum dispersion are neglected. A significant development 
provided in this paper is the dramatic simplification incurred by algebraically performing the summation over the 
spin states that are incorporated in the electron propagators. The resulting expressions in Section [I A (first stated in 
[pT|) are relatively compact, and of an appearance familiar to Landau representation/S-matrix theory applications to 
magnetized environments (i.e. including associated Laguerre functions). Furthermore, h ere t he integrations over the 
momentum parallel to the field are performed analytically for the first time in Section [IB , rendering the splitting 
rates in most amenable forms (see Eqs. [ p^ and [p^) that are optimal for numerical applications: the analytic forms 
presented consist of just triple summations over Landau level quantum numbers of the intermediate states. These 
general results are valid both below and above pair threshold at non-resonant photon energies, and provide substantial 
advances over the work of j^]; they are much more suitable for numerical evaluation since many cancellations have 
been eliminated algebraically. 

Two specializations are discussed in Section III, primarily to (partially) demonstrate equivalence of the Landau 
representation formalism presented here with extant proper-time/effective Lagrangian limiting forms for splitting 
rates, and simultaneously to serve as a check on the mathematical manipulations of this paper. Results are presented 
for all three polariza tion m odes permitted by CP invariance in the limit of zero dispersion. The first asymptotic 
regime is (see Section III A) for highly supercritical fields, B ^ Be = m^c' /eh , where in the case of -L^|| || , the limit 
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was found to concur with a recent analytic result that was obtain ed by Baier et al. ||^, while new results were obtained 
for the other two modes. In the second specialization, in Section [II B , asymptotic results for energies u) <C ttk? well 
below pair creation threshold were obtained, reproducing the cubic energy dependence of the amplitudes obtained 
by other QED techniques. Moreover, new and compact expressions for the scattering amplitudes in this low energy 
limit arc derived in terms of the logarithm of the Gamma function, its integral and their derivatives. These simplified 
forms in Eqs. (^) and (|4^) are also produced from extant integral forms for splitting matrix elements derived first in 
, thereby facilitating the first analytic demonstration of the equivalence of splitting rates obtained by the S-matrix 
formulation in the Landau representation and those derived using Schwinger-type techniques. 



II. THE GENERAL S-MATRIX FORMALISM 



The rates for photon splitting within an S-matrix formulation can be developed using a variety of conventions; 
here the Landau representation used by Mentzel, Berg and Wunner j|] is adopted, and formal developments lead to 
an independent confirmation of their analytic derivation. Specifically, for a field B — (0, 0, i?) , this approach uses a 
representation of the electron/positron wavefunctions as eigenstates of the magnetic moment (or spin) operator /i^ 
(with /X = ma + 75/9ct x [p + eA(x)] ) in Cartesian coordinates within the confines of the Landau gauge A(x) = 
(0, -Bx, 0) . Such states turn out to be very convenient because they generate useful symmetry properties; they were 
identified by Sokolov and Tcrnov p9|, who dubbed them states of "transverse polarization." Let e, e' and e" denote 
the polarizations of the initial and final (primed) photons (e, e', e" =_L, || ), and /c^ = (^i k) , k'^ — (w',k') and 
k"^ ~ (w",k") denote the absorbed and produced photon four momenta. Then the total rate for splitting via the 
polarization mode e e'e" can be written, using Eqs. (27)-(29) of [||, in terms of the S-matrix element s'^J} , which 

(31 

is the sum of six terms S jl ^ corresponding to the six viable time-ordering possibilities: 
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(1) 



where V and T denote the volume and time associated with the interaction calculation and the factor of 1/2 out 
the front avoids double counting of the final states. The priming convention adopted throughout the paper is one and 
two primes for the produced photons and no prime for the initial photon. Since the S-matrix element contains a delta 
function (5^(fc^ ^ ^ ^^i') prescribing four-momentum conservation for splitting, it is squared in the usual way using 
~^'ti~ ~^ I (2''')^] <^'*(fc/i ~^'fj,~ ^Ai') ■ -l^ote that the S-matrix element should possess a cubic dependence 
on photon energies when well below pair creation threshold, due to parity symmetry jDhoton gauge invariance, and 
the antisymmetric nature of the electromagnetic field tensor; details are discussed in 

Before writing down expressions for the S-matrix element terms, it is appropriate to identify the dimensionless 
convention that shall be adopted throughout this paper. Since the electron rest mass m is the only mass that enters 
into this QED problem, we opt to scale all energies by mc^ and momenta by mc unless otherwise specified. This 
includes a scaling of mc^/h for photon frequencies ui . In the spirit of this convention, we choose to use the symbol 
e to represent dimensionless electron energies and reserve E ( = emc^ ) to denote "dimensional" energies as in . 
In addition, the magnetic field will be expressed in terms of the quantum critical field Be = m^c'^/eh hereafter, so 
that B = 1 denotes a field of 4.413 x 10^^ Gauss. 

The convention for polarizations is identical to that assumed in 0| , who opted for real polarization vectors with zero 
time components. The polarization states _L and || are defined according to whether the photon's electric vector lies 
either perpendicular or parallel (respectively) to t he pla ne containing the photon's momentum k and the (uniform) 
magnetic field B vectors, the convention of ||l|,[l3|.|G|, |r()| . [lO[[3(i[| . In the limit of zero dispersion, three polarization modes 
are permitted by charge/parity (CP) invariance in QED, namely -L— >|||| , _L^_L_L and ||^-L|| . However, Adler Q| 
showed (see also pl| ) that for weak vacuum dispersion (roughly delineated by B ^ 1), where the refractive indices 
for the polarization states are very close to unity, energy and momentum could simultaneously be conserved only 
for the splitting mode _L— >|||| . This kinematic selection rule applies to gamma-ray pulsar magnetospheres where 
plasma dispersion is negligible. In magnetar models of soft gamma repeaters, where supercritical fields are employed, 
strong vacuum dispersion arises. In such a regime, it is not clear whether Adler's selection rules still endure, since his 
linear dispersion analysis omits higher order (quadratic) contributions p^ , p^ to the vacuum polarization tensor (e.g. 
those that couple to photon absorption via splitting) that may become significant in supercritical fields. Furthermore, 
plasma dispersion effects, which can nullify the vacuum selection rules, may be quite pertinent to soft gamma 
repeater magnetospheres, rendering them distinctly different from those of conventional pulsars. Therefore, in the 
interests of generality, consideration of all three CP-permitted splitting modes is adopted throughout this paper. 
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The derivation of the S-matrix element proceeds along lines identical to those in Mentzel, Berg & Wunner with 
the result being an exact reproduction of their analytic formalism, as reported in Weise, Baring & Melrose ; for 

details, one is referred to Of the six s'^^-'j contributions to Eq. (Q), it is sufficient to explicitly present just one: 



S 



(3) 



^ ^ nn'n" a a' a" 



P+-+(k'OP+-4kOp-+ ,(k) - p,--,(k-)p+- (kOi?-+ ,(k) 

ee'e" (e + e" + cj' - le) {e' + e" + lo - ie') 



(2) 



where a{ = e'^/{hc) is the fine structure constant, and the energies e and Eq are defined in Eq. below, with 
similar definitions for the primed energies involving primed quantum numbers and momenta of the virtual electrons. 
Using the J notation in Eq. below, 

^n'n(k") = J{-K I - 1, ~K I 71, 0) [k[*K4 + k'^ K2] c'L 

+ J{-k'; I n', -k'^ I n - 1, 0) + k'*K3] e'l 

- { J(-fc;' I n', -fc;' I n, 0) W*n, + k'^k^] - J(-fc^' | n' - f , -k'^ \n-l, 0) [k'* + 4*^i]}e'; , 
KnA^) = Ji-K I " - 1, I n", k'y) [k14* + kIkT] e'_ 

+ J{-k'^ I n, I n" - 1, k'y) [kIk'I* + <4*] eV (3) 

- {^J{-K n", k'y) [K2K2* + i^li^T] - Ji-K I »^ - 1, I n" - 1, k'y) [kIk'I* + KaKglje'^ , 
T^n'^n'i^) = Jik. I n" -l,k'y\ n', ~k'^) [k'^k'^ + 4<] e_ 

+ J(fc, I n", k'y\n'~ 1, -fc;') [4< + 44'] e+ 

- { J(fc, I n", k'y I n', [44 + 44] - J(fc, I n" - 1, 4 | n' - 1, -k'^) [44 + 44]}e, , 

where the polarization vector = {0, Cx, Sy, ez) is specified by e± — Cx i isy and , and similarly for the final 
photon polarizations (primed) . The other three I? s in Eq. are not displayed here for brevity; they can be obtained 
from those in Eq. (|3|) simply by the interchange e+ <-> e_ of polarization components (and similarly for primed 

components) and a relabeUing of the Js that produces a correspondence I?^,^(k) T^nn~'' (k) • 

Several notations need to be identified. First, the particles have energies e , and momentum components Pz along 
the field. The energies e and eo that appear here are, respectively, with and without the parallel momentum pz '■ 



e ^ ^/l+pl + 2nB 



1 + 2nB 



(4) 



with n denoting the Landau level quantum numbers, as usual. The other quantum number pertaining to the 
eigenstates of fiz is a — ±1 , which signifies the spin state of the fermions ( [Q used the label t ; here the notation of 
Melrose and Parle |f5| is preferred), and satisfies Hztp = aeoip . It does not appear explicitly in e, but is embedded 
in the spinor coefficients Ki : 



K2 



/ s_ s+ 

S+ 










-6+ 



\ 





( 4(£o + l)(e + eo) \ 



£0-1 



£0 + 1 



(5) 



for i5+ = 5a,i and i5_ = (I — 5nfi) 5a,-i where the spin quantum number a takes on two values except for the n = Q 
ground state (zeroth Landau level), where only cr = 1 is permissible. Here Sij is the familiar Kronecker delta. The 
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primed coefBcients k[ and k" are similarly defined in terms of primed momenta and Landau level quantum numbers, 
subject to the momentum conservation implicit in Eq. (||). 

The J functions that appear in Eq. (|^) are integrals over the oscillator functions (Hermite polynomial products), 
a form undeveloped in Here, Eq. (7.377) of is employed to express these integrals analytically in terms of 
generaHzed Laguerre polynomials, U^~^{x) (see also Eq. (47) of |^): 



J(a|n', /3) = exp( ^ [/? + /3']) {zg-*^}" (6) 

where p ( > ) and the phase ij) are introduced for convenience of notation: a — p cos "0 and j3' — (3 = p sin ij) . 
Note that the /3s are always kyS and a is always a . Here the In',n functions follow the Sokolov and Ternov 
convention |29) up to a factor of nl , being related to the J functions of Melrose and Parle |l5j, and both are defined 
in terms of the generalized Laguerre polynomials (see js^): 



In'A^) = (-!)"'""/«,«' (:r) = K'-ni^) ^ \ ^ e-^/2a;("'-")/^Lr;-"(x) , n' > n . (7) 




Values for n > n' are obtained by interchanging indices, as indicated. Hereafter, the (modified) Sokolov and Ternov 
convention for writing the Laguerre polynomials will be adopted. Complex conjugation of Eq. (^) can be used to 
establish the identity 



n'n 



(k) = (-i)"'-"{i?-«r«'(k)}* , (8) 



noting that the k products are either purely imaginary or real, for all choices of spin quantum numbers. The three 
factors like (—1)" ~" appearing in the second product of three I?s in Eq. (||) cancel, leading to this product being 
just the complex conjugate of the first three Vs. This useful symmetry property clearly underlines the convenience 
of the Sokolov and Ternov choice of wavefunctions when adopting real components for the photon polarization. 
The form of the contribution to the S-matrix element in Eq. (^ is identical to that for S^j^^\ given in Eq. (25) 

of Mentzel, Berg and Wunncr 0. In the same fashion, it can be found that the expression derived here for Sjl\ 
is absolutely identical to Eq. (26) of [Q, thereby providing confirmation of their analytic developments; it can be 
obtained by using the substitutions k' ^ k" , n' ^ n" and e' ^ e" in Eq. (|). All other sfl^ contribut ions result 
from application of the cyclic permutations 

p 7// // h' ' h' ' h * h * h" " ■ 



P ^ ■ ~h p* —y h' p' h' p' ^ h" p" h" p" ~h p* 



(9) 



where A;^ = (cj, k) , = (0, Cx, Sy, Cz) , etc. Observe also that a minus sign and the complex conjugation of the po- 
larizations are always associated with the initial photon since it is absorbed in the process. Given these permutations, 
the crossing symmetry for splitting is manifested in the following relationship between the various terms like those in 
Eq. (|) that contribute to Eq. (|l]): 

c(3) _ p q{3) cf(3) p q(3) q(3) _ p q{3) q{3) p ri(3) /in"! 

where the permutations act as operators. This symmetry can be expressed in a multitude of ways using the identities 
P+iP_i = / = P_iP+i and P|i = / . 

It is important to remark that the derivation of analytic forms by Mentzel, Berg and Wunner is not the first in the 
literature relating to S-matrix applications to photon splitting. The papers by Melrose and Parle dealing 
with various aspects of QED in strong magnetic fields, specifically from a wave dispersion/response tensor approach, 
constructed the S-matrix element for splitting in Eqs. (46) and (47) of |p^ , which incorporated the quadratic vacuum 



response tensor given in Eq. (36) of |14|. This tensor is obviously of a standard S-matrix Landau representation 
appearance. Eqs. (^) and ^ can be generated directly (and also S^^^2 ) from the Melrose and Parle evaluation after 
a modicum of algebra. Hence, Eqs. and here, and Eqs. (25) and (26) can be used as reliable starting points 
for further S-matrix developments. 
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A. Analytic Reduction: Summation over Spin States 



The form in Eqs. (|^) and (^) is quite cumbersome. It can be simplified considerably by (i) specializing to specific 
but representative directions of photon propagation and (ii) analytically performing the summations over spin states 
a , a' and a" . Restricting the photon motion to the x -direction yields photon motion perpendicular to the field: 
since splitting is coUinear in the non-dispersive limit discussed earlier in this paper, it follows that fcz = fc^ = fc" = . 
This choice dramatically simplifies coefficients of the Laguerre polynomials in Eq. (|3|). Without significant loss of 
generality, setting ky — k'y = ky — removes nearly all of the phase factors in the definition of the Js in Eq. (|^), 

leaving just z" ~" . Three such factors emerge in the triple product of Ps, leading to a factor of (—1)" ~" . 

The CP symmetry possessed by the splitting process becomes most evident at this point, since it is now simple to 



derive the CP selection rules. The specification of ky = k'y = ky = and = k'^ = k'!, = yields only one possible 
component of polarization perpendicular to the field, = Sy — ~ie+ — ie- and one conceivable component of 
polarization parallel to the field, ey = Sz (and similarly for primed quantities). The polarization (electric field) 
vector of the photons is, of course, normal to the photon momentum vector, which automatically spawns the notation 
for the two possible polarization states: _L : e± — 1, = and || : ej_ — 0, e|| = 1 . From the presence 
of subtractions in the numerators of the integrands of Eq. ^) together with the complex conjugation property in 
Eq. (R|) and the proportionality of the I? s to factors like ^" , it follows that only terms with an odd number of 



e± 



("3") 

factors contribute to Sj:^\ . 



i.e. terms proportional to en eye 



and 



e±e', e'\ 



All other terms 



(3) 

this is also true for all other S 



It 



cancel identically to zero. By virtue of the permutation symmetries in Eq. (|^), 

is then trivial to deduce the CP selection rules for photon splitting, namely that the only permitted transitions are 

, , 11-^11 • (11) 

The three other splitting transitions all have S-matrix elements that are exactly zero for collinear photon momenta, 
and hence are forbidden. This technique for CP selection rule derivation was implemented in These restrictions 
are simply consequences of the charge conjugation (C) and parity (P) symmetries of the splitting process, i.e. relating 
to the transformations k — k and B — B . 

The summation over the spin states a , a' and a" ( = ±1 ) produces a dramatic simplification in the appearance 
of the S-matrix elements. Such spin summations act only on the products of the s that appear in Eq. (^; the 

. . (3) 

algebra is lengthy but straightforward, being facilitated by pairing Sj:/^ terms with denominators that differ only in 
the sign of their photon energies. The total splitting rate in Eq. (ol) can be written in the form 



Re — >e 



2^ 



mc 
IT 





Me^e'e" 







(12) 



where 



is implicitly understood from the conservation of four-momentum. While these rates will be 



expressed for photon propagation normal to the uniform magnetic field, the results for general photon obliquities 



cj smt 



together with 



to B can be obtained via a simple Lorentz transformation: uj — *■ uj sinO , u 
an extra multiplicative factor of sin 9 applied to the rate in Eq. dlj) . 

The momentum dependence in the integrands of Eq. (|^) can be simplified by forming sums of the products of energy 
denominators. Separating such sums into real and imaginary parts via the representation — ^ ; leads 

to the definition 



yR _ 



+ A 



-77)V 
1 



1 



' to) 



1 



(e' + e" - cj) (e + e' - uj") ^ (e' + e" + u;) {e + e' + u;") 
ie + e' + Lj") (e" + e ^ u') ^ [e + e' ^ u;") (e" + e + u;) 



(13) 



of the real part, where A and fj, assume the values ±1 . The imaginary part is not explicitly stated since it will not 
be of use in the subsequent developments. The momentum integrations over these then assume one of the forms 



T — 



Pz dPz yK 



(14) 



J — I ^1,-1 ' ~ / ~p~ -1,1 ' ~ / ^-1,-1 
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for n = or 1 ; generalizations to complex S^.^ (relevant to calculating splitting rates above pair threshold and near 
pair resonances) are routine. These manipulations yield the following compact forms for the A^e^e'e" coefficients in 
Eq. (P): 



M 



i^llll 



V2n"B 



J"+Ii-Io 



A 



■nil 



A 



i^llll 



M 



M 



B_ 

4 



4 



^ (-l)""-"'|V8nn'n"B3Xo A^^ 



\/2nB 



A 



i^llll 



V2n"B 



J" — 2^1 — 2o 



A 



-\/2n'B 



A 



(15) 



^ (-l)""-"'<{ \/8nn'n"B3Xo Af^^" +\/2n"S 



^"+Xi-Xo 



A 



IHJ-II 



A 



IH^II 



V2nB 



J-Ii 



A 



IH^II 



results that are to be used in conjunction with Eq. (|l^). The factor of — _B/4 is introduced to render the scaled 
amplitudes positive, and also to afford a direct mapping onto limiting forms obtained P,fO[ by the proper-time 



technique, as will become evident in Section [II. The A,^^*^ are differences of triple products of generalized 
Laguerre polynomials (defined in Eq. 0); for _L^||| 



(16) 



A--" 


= /" 


l,n' 


Jl 

-1 ^n' 


,n In 


.n" — 


T" 


, /' 

n' n 


'-l,ri 


-1 In'-l,n" 


A^" 


= /" 


l,n' 


Jl 

-1 ^11' 


-l,n- 






-I" 


,/', 

n n 


In' ,n"-l 




n- 


l,n' 




-l,ri- 


-1 In 


,ri"-l 


n 


,1', 

n n 


In'-l,n" 


Ar" 


= /" 


\.n' 


J' 

-1 ^n' 


,n In 






n' n 


'-l,n 


— 1 In' .n" — 1 



where the Sokolov and Ternov representation of the associated Laguerre functions in Eq. 
the priming notation 



In' ,n 



I ' I — 



I' , 

n' ,n 



In' ,n 



.'l2 



2B 



I", = I , 



,//l2 



2B 



is used together with 
(17) 



thereby aiding brevity. For the _L^_L_L mode. 





n,n' 


T' 

-1 ^n' 


,n-l In' ,n"-l 


-I" 


J' 

l,n' n' 


-l.n In' 


-l.n" 


^2 


= /" , 


-ll'n' 


-l,n In'-l,n" 


-I" 

71- 


J' 

l,n' ^n' 
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(18) 



and the results for the ||— >-L|| mode are not explicitly stated since they can be obtained by exploiting crossing 
symmetries: the inverse of the permutation in Eq. ([l^) yields the transformation Aj*"^"" — Af^^" , A^^"" — > 

The A^^*^'*^" can alternatively be expressed using the 



A 



IH^II 



, A 



-L-llll 
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IH^II 



, A 



-L-llll 



A 



IH-LII 



functions of Melrose and Parle as in [|ll|. Note that the potential subtlety of having to include factors of 1/2 for 
some contributions from ground intermediate states is eliminated by the specific choice of the Sokolov and Ternov 
wavefunctions. 

The comparative simplicity of the reduced form of the S-matrix element relative to Eq. (^ is both notable and 
comforting. Unlike Eqs. (25) and (26) of ||l|], this developed form of the splitting S-matrix element has an appearance 
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familiar to S-matrix applications of QED in the Landau representation to strongly-magnetized systems, with products 
of generalized Laguerre polynomials multiplied by simple combinations of energies and momentum components. Ex- 
amples of previous work bearing such familiar forms focus la rgely on lower-order QED processes and include studies 
of synchrotron radiation [ ^J30| , single photon pair creation , and vacuum [p4| and plasma [^5| polarization. 

For the purposes of the analysis in the next section, it is pertinent to define the cyclic permutations 

. , „ , (19) 

n n , n— >n, n — >n, 

in the spirit of the P+i permutation in Eq. (^. These permutations will appear repeatedly in the developments 
below, and lead to the following transformation properties of Eq. ( p^ ) : 

with "Sj^i being invariant, symmetries that are consequences of the arrangements of electron and positron propagators 
in the Feynman diagram for splitting. These translate into obvious mappings between J , J' and J" and an 
invariance of the /„ . It is also easily seen that under this cyclic permutation, the factor in braces in the summation 
for is invariant, while the equivalent factor in the summation for maps over (up to a minus sign) 



to the factor in braces in the A^^^im summation. As will become evident in Section III, the remaining powers of 



- 1 in the summations do not provide any unsatisfactory interference in the limits of low photon energy ( w ^ 1 ) 
and high fields {B ^ 1 ) , so that permutation symmetry can be extended to the total amplitudes in these specific 
parameter regimes. 



B. Analytic Reduction: Integration over Parallel Momentum 

Further analytic development is not only possible, but also desirable, given that the integrations over the momentum 
Pz parallel to the field can be expressed compactly in terms of elementary functions. Such tractability facilitates both 
numerical evaluations and the derivation of asymptotic limits. In proceeding, since results are sought at energies 
sufficiently remote from pair creation resonances, the imaginary parts of the denominators in the Sa,/j are dropped 
in all further considerations, i.e., we consider only the functions = ReE^,^ . 

It turns out that carefully-constructed contour integrations in the complex pz plane do not facilitate the pz 
integrations. Hence the first step in integrating over pz is effected by the more cumbersome and less elegant approach 
of completing the squares and rationalizing the denominators using products of factors like (e' ± e" ± lo) . These 
factors define poles pij of the Pz integration for i and j being some combination oi n , n! and n" . Such poles fall 
into two types: pair creation ones (e.g. see 0) that contribute only above pair threshold, due to the structure of the 
splitting rate, and cyclotronic ones that must be considered below pair threshold. The appearance of such cyclotronic 
poles is an artifact of the rationalization of denominators, so that they are really pseudo-poles of the subsequent 
analysis; a consistency check on the algebra is that the S-matrix element be effectively continuous across them. It is 
convenient to define energies that correspond to the pij poles: 



{J'f+M-M' _ <J+M' -M" _ {uj'f+N"-N 



(21) 



and three others paired with these, which are obtained via the relations e„'„ -I- £„„' = w" , £„"«' -t- £n'n" = and 
£nn" + £n"n = ^' ■ Here the notation 

A/" = 1 + 2nB , N' = 1 + 2n'B , N" = I + 2n"B (22) 

is used for the purposes of abbreviation. Observe that, taking advantage of the subjectivity of such definitions, a 
minus sign appears in front of the expression for en'n" , a choice that preserves symmetries induced by the mapping 
in Eq. ( |l9| ) in the results that follow. These definitions spawn the following useful identities for the momentum poles: 

■ - M" (23) 

-M , 

which immediately imply the possibility of poles along the imaginary axis. In fact, > — minjA/", Af'} , with 

equality for uj" = \H — , and likewise for the other poles. Note that for the one- vertex calculations of cyclotron 
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emission and single photon pair creation and annihilation, the requirement that such poles be reaLcorresponding to 
real components of particle momenta on external lines, is precisely what generates thresholds (e.g. [ |l7[ ) and kinematic 
cutoffs (e.g. [Q) for transitions involving various states. 

The rationalization of the denominators yields relatively compact decompositions for these sums, after much can- 
cellation and simplification. They take the form 



the simplicity of which is contingent upon the energy-conservation restriction uj" — uj ~ uj' . Here 

W = ujJuj" + -uj'N' -uj"N' . 



(24) 



(25) 



Identities such as = —2ujuj'{en"n' + £n"n) prove useful in the ensuing analysis. The c\^^ and t\^^ coefficients 
assume simple forms when expressed as partial fractions. Consider first the result for , which has the coefficients 
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(26) 



Observe that a cyclic symmetry is immediately apparent: 'S^i is invariant under the permutation in Eq. (p^), as is 
evident from its original definition in Eq. (p^). Similarly, the algebraic developments yield coefficients for the _i 
sum, which appears in the ee'A2 terms, as 



c-i.-i = 
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UJLU 



^n"n- 
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(27) 
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^n"n 
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The coefficients for the sum TPi^ that appears in the ee" A^ terms and the coefficients for the sum that 
appears in the e'e"A4 terms are similar: there is little need to state them explicitly, since the coefficients possess a 
relationship to each other due to the permutation symmetry enunciated in Eq. ([2C|). 

Given these decompositions, it is now fairly straightforward to evaluate the integrations over , expressing them 
in terms of the an elementary function / with real arguments Sij : 



f{Af, £) = P 



dpz 



£ 



1 



./AfT^_ - M - Pi 



= < 



^£-^-N 
2 



^Af-£^ 



arctan 



£-^£^-N 



if 52 > , 
if < £2 < ^ 



(28) 



for real £ . The identity arctanz = (l/2i) logg[(l-|-«z)/(l — zz)] with z = —£I\JM — f ^ has been used to map across 
the singularities at £ = ±y/Af (cyclotronic below pair threshold) and guarantee bounded and continuous behaviour 
of f{Af, £)/£ at £ ^ . The integral identity in Eq. ( p8| ) can be established quickly with the aid of result 3.513.2 in 
p3| , using the substitution pz = VAT sinhi and partial fractions. Note that real values (either positive or negative) 
of £ are guaranteed by the formalism here, with £ — being improbable due to the discreteness of the quantum 
numbers n , n' and n" . 

The integration of the coefficients Iq of the Ai terms for each of the polarization modes are then straightforward, 
and the identities in Eq. (E3h can be used to advantage. Similar terms appear in the Xi integrations of parts of the 
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coefficients of the other Ai terms, which also possess integrands with terms proportional to 1/e, 1/e' and that 
formally lead to divergences that cancel each other (an artifice introduced by the rationalization of the denominators). 
Using partial fractions, the divergent contributions can be written as integrals over the finite range — p < Pz < P , 
rearranging to subtract off exactly-cancelling terms, and then taking the limit as p — > cx) . Similar manipulations 
are used for the J" integration over E^j^ _i , where again the leading order terms are individually divergent yet 
collectively convergent. Partial fractions can again be used to enable rearrangements and separate the divergent 
terms, which are then integrated over finite ranges as with the Xi evaluation. The results are encapsulated in the 
identities 

Xi = £■ + jjT < J-nn' + Pn'n" n'n" + Vn"n ■f'n"n \ 1 (29) 



where 

Tnn' = /(AA, Enn-) + /(A/"', En'™), ^n'n" = /(A^', £n-n") + /(A^", ^n-n = £„"„) -f /(A^, £„„"), 

(30) 

and 

4-jT logeAA- -i- log,AA' - log,AA" . (31) 
No further integration is necessary: the cyclic permutations in Eq. ( p9| ) can be used to quickly derive expressions for 



J' and J from Eq. (29) 



At this point, it is salient to remark that the divergences at E"^ = N in the functions f{J\f,£) pose no problem 
for the integral evaluations in Eqs. (p9|), because these functions always appear two at a time. Below the pair 
threshold, these divergences are cyclotronic in nature, being encountered when u) \VM' — VA/""] or for similar 
circumstances for the other photon energies. As w tends to such a limit, for example, we observe that ^/ Af' 

and £„"„' —\/Af" when Af' > N" (without loss of generality). This opposition of signs guarantees cancellation of 
divergences when the arctan form of f{J\f,£) is used (arctan(l/z) Tr/2 — z as z ^ 0), so that continuity across 
cyclotron "pseudo-resonances" emerges naturally from Eq. ([29|), consistent with the continuity of the functions. 
Continuity across pair resonances does not arise above pair threshold, so that true divergences emerge. 

The incorporation of Eq. ( p9| ) into the scaled matrix elements in Eq. ( p^ constitutes the final product of the general 
analytic developments in this paper, providing rates valid for all energies below pair threshold (and applicable for 
non-resonant energies above threshold), and for photon propagation normal to the uniform magnetic field. They are 
eminently suitable for numerical computations, having improved upon the analytic formalism of Mentzel, Berg and 
Wunner Q| (i.e. Eq. j2j) by performing the summations of the spin states and integration over the momenta parallel 
to the field that are associated with the electron propagators. Such developments are prudent prior to numerical 
evaluations due to the large degree of cancellation in these sums and integrations. 



III. ASYMPTOTIC LIMITS FOR HIGH B OR SMALL a; 



A fruitful extension of this analysis is the exploration of the simplification of the scattering amplitudes and rates in 
two particular asymptotic regimes, namely the limit of highly supercritical fields, i? 3> 1 , and the specialization to 
photon energies well below threshold, i.e. w ^ 1 . The benefits of such an investigation are twofold. First, it provides 
the first unequivocal analytic demonstration of the equivalence of splitting results from the S-matrix formulation in the 
Landau representation and effective Lagrangian/proper-time results from Schwinger-type formalisms in well-defined 
parameter regimes. In doing so, it serves as a powerful check on the developments here. Second, in the w <C 1 case, 
it identifies a new, satisfyingly compact representation of the scattering amplitudes in terms of special functions that 
leads to an efficient means of computation. 

These two parameter regimes are encompassed under the single limit cj^ ^ 1 -I- IB , which thereby identifies 
the appropriate series expansion of the generalized Laguerre polynomials that appear in the amplitudes. For small 
arguments x , the leading order terms in the series for In',n{x) can be found in the Appendix of ||l5[. Given that n , 
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n' and n" cluster in a manner such that |n' — n| ~ \n" — n\ ~ 1 , this series converges rapidly provided nx ^ \ . 
Hence nuj'^ /{2B) actually represents the true expansion parameter here, with to' and oj" being similarly bounded. 
The leading order terms of such expansions for the Af"**^ are linear in the photon energies, while the next higher 
order terms are cubic; a more detailed exposition can be found in Weise, Baring and Melrose |]ll|. The series for the 
integrations of pz , namely , Ii , J^" , and (which do not depend on the polarization mode) are expansions 
in cj^/(l + 2B) rather than u)'^ /{2B) . They are independent of photon energy to leading order, with a quadratic 
scaling with energy to next order. The series for <C 1 + 2B possess logarithmic character in the quantum numbers 
in situations when no two of them are equal (i.e. M ^ M' ^ M" ^ M ): 



Xi « -J" 



41og,AA , 41og,AA^ , 41og^AA" 

{M - M'W -M)^ {N' - N"){N - N') ~^ {N" - - M") ' 

2N\og^N 2N'\og^N' 2N" log^N" ^^^^ 

{N~N'){N" -N) {N' - N"){N - N') {N" - N){N' - N") ' 



and additionally involve inverse trigonometric functions when two n s (e.g. for M = M' ) arc in fact equal: 

(33) 



4 , N 4 



r 2M" M 2 , AN-{Lo"f iu" \ 4 uj" \ 



{Af- 



where 



Q(a;) = arcsmx ^^^^ 



z\f\ 



and the identity arcsina; = arctan[a;/Vl — has been invoked. This retention of the inverse trigonometric functions 
is particularly relevant for determining the high B limiting forms of the scattering amplitudes. Relations similar to 



Eq. (33) exist for M = M' and M' = Af" , obtained by the cyclic permutations through Afs and photon energies. 
The lengthier higher order (quadratic) terms are not explicitly stated for the sake of brevity. This concludes the 
preamble that guides the reader in the subsequent specializations. 

A. The Special Case of B > 1 

This regime is of particular relevance to the study of magnetars such as soft gamma repeaters. For the two modes 
-L— >|||| and ||— >-L|| , only the leading order terms for the Ai and the momentum integrals presented in Eqs. ( |3^ ) 
and (p3) are required. Consider first the reduction of A^^^im . Here the A2 and A3 terms contribute leading order 
terms only through n" = 1 , n — n' = and ri' = 1 , n = n" ~ cases, respectively, where it is necessary to use the 
full forms in Equation (^), and inverse trigonometric functions appear through the Q{x) function, which assumes 
the arguments x = uj' /2 and x = u" /2 . A similar n = 1 , n' = n" = term is identically equal to zero by virtue of 
the A4 factor. The contributions from the Ai and A4 terms possess an entirely different character, being infinite 
summations over n , with the values of n' and n" being constrained by |n' — nj + |n" — n| < 1 , producing five 
groupings of the indices. The series is evaluated by truncating the sum at n < k , relabelling one of the logarithmic 
terms, and then taking the limit k 00 . The net result is (for ui < 2) 

A^.^llll - arcsm(^j + ;jr^^=^ arcsm(^ j - , S » 1 , (35) 

which, when combined with Equation (|l^), yields the asymptotic high-B result derived by Baier et al. and 
reproduced independently by Baring & Harding [|l^; the overall rate for -L— >|||| approaches a value independent of 
B . Observe that the manifestations of the pair creation threshold for each of the final photons of || polarization 
(i.e. at w' = 2 and uj" — 2) are the individually-divergent coefficients of the inverse trigonometric functions. Yet, 
collectively, due to the energy conservation relation lo = u' + lu" , such divergences cancel each other to yield a finite 
overall result as a; ^ 2 . For the incident photon of _L polarization, the pair threshold of 1 + Vl + 2B is remote from 
w = 2 so that it would only become explicitly apparent when the amplitude was evaluated to higher order in B . 
Note also that the w <C 1 limit of Eq. (^) is ujlu'lu" /6 and reproduces results obtained in Q and The functional 



form of Eq. (pSD is plotted in Figure 
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scaled by co^ , on the fractional energy lo' /uj of one of 
(in units of mc^ ), as labelled. Only the two polarization 



FIG. 1. The dependence of the scattering amplitudes for B S> 
the produced photons, for three different incident photon energies lj 
modes with amplitudes asymptotically independent of B (in units of Be ) in this ultra-quantum limit are depicted, namely (a) 
_L— >|| tl and (b) [|— >-L [| ; their functional forms are given in Eqs. ( ^ ) and (^^, respectively. The shape of the amplitude curves 
for _L— >±± is independent of u) and is very close to that of the uj = 0.1 curves in panels (a) and (b). While the ^^|| || curves 
are necessarily symmetric about ui' — uj/2 , asymmetry is present in the [|— i-J- || case where cu' represents the final photon of 
_L polarization. Note that the magnitude of M\\^± \\ diverges as pair threshold w = 2 is approached. 

The equivalent result for the splitting mode requires little additional algebra given that it can be obtained 

from the analysis just above using the cyclic symmetry transformations of Eq. (|l^). Carefully keeping track of signs 
and all photon frequencies by relabelling at the beginning of the manipulations, the roles of the A4 and A3 terms 
are interchanged, and the obvious result emerges: 



M 



IH^II 



arcsm 



4a; 



arcsm I 



B > 1 



(36) 



While not established before in the literature, the low energy limit of this, namely « ojuj'uj"/6, yields the 

differential rate from previous expositions of low energy approximations. The form of Eq. ( p6| ) is displayed in 
Figure ^ exhibiting the asymmetry expected under interchanges uj' ^ uj" . In this case, pair threshold structure in the 
amplitude appears again for the two photons of parallel polarization (i.e. at uj = 2 and uj" = 2), and is also absent 
for the produced _L photon, being of higher order in B . Consequently, the amplitude possesses a real divergence at 
w = 2 , a noteworthy occurrence that is illustrated by comparing the two panels of Figure ^ Such divergences, which 
are not integrable over uj (and therefore patently different in nature from the resonances encountered in rates for 
7 — > ), are characteristic of the photon splitting rate near resonances at and above the pair threshold of uj — 2 , 
corresponding to the creation of virtual pairs in various excited states. In fact, near such resonances, photon splitting 
necessarily becomes first order in at like pair creation as the intermediate states "go on-shell." 

The rapid increase of the rate of ||-^-L|| relative to that of -L^|||| is exhibited in Figure |[ where the rates have 
been scaled by the low energy (uj <^ I) limiting forms {R{uj) oc uj^) discussed in the next subsection. This particular 
scaling is chosen to illustrate deviations from the uj ^ 1 asymptotic forms, and therefore to demonstrate the need 
for relinquishing use of them when samplingphoton energies near pair threshold, a parameter regime very relevant to 
certain astrophysical calculations (e.g. see p7| , p8[ |). The dominance of the i?||_|._L|| over i?_L^|||| near uj = 2 apparent 
in these B ^ 1 results becomes substantive in parameter regimes where the weakly-dispersive vacuum (i.e. for 
B ^ I) polarization selection rules for splitting derived by Adler [Q (which prohibit ||— >-L|| and _L^_L_L splittings) 
may not apply if non-linear contributions to vacuum polarization or plasma effects are significant. This underlines 
the saliency of a detailed determination of the dispersive properties of the magnetized vacuum or plasma medium 
appropriate to a particular astrophysical scenario. 

The derivation of the B ^ 1 form for the amplitude for _L^_L_L differs significantly from the results just expounded. 
First, contributions from n" = I , n — n' — and n' = 1 , n ~ n" — Q and n = 1 , n' ~ n" = Q combinations are 
identically equal to zero by virtue of each of the associated factors. This automatically implies that no inverse 
trigonometric functions that have arguments independent of B appear in the amplitude, a property not possessed by 
the other splitting modes. The consequences of this are twofold. First, this cancellation implies that the scattering for 
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_L— »_L_L is of a higher order in B than for the other two sphtting modes. Second, since any potential appearance of 
inverse trigonometric functions spawned by the forms in Eq. (^3|) involves arguments that depend on B through the 
Ms, these arguments are always small when B ^ 1 , precipitating a redundancy with the low energy limit. Hence, it 

follows immediately that the scattering amplitude for _L— »J L in the regime of highly super-critical fields is identical 

to that of the B ^ 1 specialization of the low energy (cj ^ 1 ) limit. As the latter has been derived in various papers 
in the literature (e.g. see [^,p|Jll|] and the subsequent section), here it is sufficient to merely state the result: 

-^^-^^ « ^ , S»l . (37) 

This extremely simple form differs profoundly from those of the other two modes because of the absence of photons 
of II polarization in the interaction. Hence any signatures of the pair threshold of 1 + y/T~\~2B of _L photons are 
absent in the domain of lj < 2 , and a scaling-type form with obvious cyclic symmetry emerges. 




.5 1 1.5 2 



0) 

FIG. 2. The total rates in the B ;3> 1 limit for the modes -L-^|| || and |[^_L |[ , computed according to Eq. using 
the amplitude formulae in Eqs. (^^ and (|36|), divided by the rates that would be computed when taking the low energy 
(a; <C 1 ) limit of these amplitudes, i.e. M_i_^\\ \\ ~ louj'lo" /Q ~ A^||_,x|| . Deviations from such low energy approximations (i.e. 
R{ijj) oc ), while significant for -L^|| || , are dramatic for ||^-L || near pair creation threshold w = 2 . 



B. Approximations for uo 



The low energy limit w ^ 1 is of interest not only because it was the regime where compact analytic expressions 
for the splitting rates were first obtained (^-IJ], but also because the analysis that follows derives simple and ele- 
gant representations of the scattering amplitudes in terms of well-known special functions that provide a convenient 
alternative option for numerical evaluations. 

The amplitudes for each of the splitting modes should exhibit a cubic energy dependence [Q when a; <C 1 . Hence, 
a necessary product of the Landau representation formalism is that terms linear in photon energies should contribute 
exactly zero. For the polarization modes _L— >|||| and ||^-L|| , whose amplitudes are identical in the low energy limit 
the demonstration of this is not dissimilar to the i? ^ 1 analysis. The w <C 1 restriction generates a single 
infinite series in n due to the clustering of n' and n" around n . The ensuing algebra in the simplification of this 
series is moderately lengthy, and requires re- indexing of the logarithmic terms to assume forms involving logg [l-|-2n_B] , 
and also some relabelling of the rational functions. Care must be taken in these rearrangements due to the infinite 
nature of the series, and the technique adopted is outlined just below. The terms linear in photon energy result in 
zero, as expected: for more details, the reader is referred to ||ll|. The next order contribution is cubic in energy, 



13 



as desired, with terms coming from a mixture of (i) the hnear terms of the combined with the quadratic higher 
order terms of the Pz integrals, and (ii) the cubic terms in conjunction with the leading order (constant) terms 
from the pz integrals. The algebra is straightforward, but lengthy and tedious, generating an exact cancellation of 
all but terms proportional to lolu'lo" . This approach leads to a reproduction of the Cj listed in entirety in Appendix 
B of Weisc, Baring & Melrose for both the -L^|||| and _L^_L_L modes of splitting. Hence there is little point in 
replicating these expressions here; the reader is referred to [0 for details. 

These results are expressed as single infinite series in the label n , which sometimes starts at n — Q , and sometimes 
begins at higher integer values (up to 3). Hence, an aesthetic goal is to rearrange some of these series so that the 
summations in each contribution begin at n = , and then add the terms in the series together. This is a non- 
trivial exercise, given the divergent nature of the series in many of the individual contributions. Hence, considerable 
care must be taken when performing the rearrangements, for which there is no unique prescription. One choice for 
relabelling the sums is adopted by Weise, Baring & Melrose though their end results expressed in their Appendix 
C do not facilitate analytic development in the most expedient manner, and were in fact erroneous (discussed briefly 
below) . An alternative and preferable choice for rearrangement of the multitude of series over the label n is adopted 
here, outlined as follows. Inspection of the various Cj contributions in Appendix B of reveals that they always 
consist of three types of terms: (i) logarithmic ones proportional to logg[l + 2{n + l)B] , for I ^ 0, ±1, ±2, ±3 , (ii) 
rational functions of 1 + 2(n + l)B for I — 0, ±1, ±2 , and (iii) polynomials in n . A unique method for rearrangement 
is to truncate all series to finite ones with n < k , and then perform relabellings so that the first two types of terms 
consist only of logg [1 + 2nB] terms and rational functions of 1 + 2nB . This approach provides no particular focus on 
series that originate with labels n > , but requires careful accounting of the remainder terms at the upper and lower 
ends of the sums, for which significant cancellation arises. The coefficients of the logarithmic functions, originally 
cubic in n , reduce to linear functions of n in this development. The consequent simplification of the series terms is 
counterbalanced by the transferral of complexity to the constant remainder terms, which are purely functions of k 
and B . Taking the limit of fc ^ oo achieves the desired (and convergent) result. 

After considerable algebra collecting together all the constituent series in Appendix B of , and performing the 
rearrangement as just prescribed, one arrives at the following series representation of the scaled scattering amplitudes: 



= UJUj'oj" lim i Te-^e'e"in, B) +TZe^e'e"{k, B) 



(38) 



for oj ^ 1 , where 



|(n, B) 



4n 3 
B^2B' 
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2B 
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2B^ ) l + 2nB 3(l + 2nS)^ 



T^^^^(n, 5) = ^ log. 



2B 



1 



1 



(39) 



2B^ l + 2nB 2B^ {1 + 2nBy 



defines the series terms. The remainders are quite lengthy, and are listed in Appendix A. 

Consider first the polarization mode _L— *JJJ| . While possibly only marginally simpler than Eq. (CI) of 11|, the 
series and remainder in Eqs. (|3^), ( ^9[ ) and (Al) naturally enable the development of a special function representation 
of the scattering amplitude. The finite summation over terms like {x + n) logg(x + n) in Eq. ( p9| ) can be expressed 
using result 44.1.2 of in terms of an integral of the logarithm of the Gamma function. At this juncture, the 
analysis begins to image parts of that generated in expressing the polarization properties of a magnetized vacuum 
via effective Lagrangian or proper-time techniques ||3^-|39[|, as should be expected. Hence, it is appropriate to adopt 
definitions from such literature as much as possible. Following [|37y38[ , here a definition for the generalized Gamma 
function ri(x) of 



log,ri(x) = ^ dt \og^T{t) + ^x{x- I) -^l0g^27T 



(40) 



= 1 , are discussed at length in pffl and outlined in 



is adopted. Properties of this function, which include ri(l) 
Appendix B. 

Using Eq. (Bl), one soon arrives at an expression for the scattering amplitude in terms of a handful of special 
(polygamma) functions, namely ri(a;) , and loggr(x) and its derivatives. This representation consists of two parts, 
one independent of k , and one that involves a limit as k oo of the remainder in Eq. ( [Al| ), combined with several 
terms incorporating the special functions with arguments that depend on k . In evaluating this limit, most terms can 
be handled in a straightforward manner, and standard asymptotic series (e.g. see [p3[) for loggr(a;) and ip{x) as 
X ^ oo prove useful. However, the treatment of the term involving the function logg ri(l + fc + 1/2B) that appears 



14 



in the limit contrib utio n is non-trivial. A series representation for this function for large arguments is required, and 
is presented in Eq. (B4). Assembling the various pieces, the limiting result as fc — > oo is 



M 



nil = ujw'uj" 



|^iogeri(2^)-^iog,r(2^) 



2BJ 



12B 



y^^ijB) -i+w-'-w+ i-o^s^^.- 1 - 3iog„2i?) 



(41) 



for a; ^ 1 . This is the sought-after compact analytic form that is comparable in simplicity to the one-loop effective 
Lagrangians calculated in [^^. Using series and asymptotic expansions for all the special functions present, it is 
routine to establish that Mj_^i\ w (265^/315) uJLu'uj" for _B <C 1 , while for B ^ 1 , one finds A^j^^im ~ ujuj'oj" /Q , 
a result obtainable from Eq. ([3^). 

The developments are similar for the _L^_L_L mode: this representation again consists of two parts, one independent 
of k , and one that involves a limit as fc — > oo of the remainder in Eq. ( |A2| ), combined with several terms incorporating 
polygamma functions with arguments that depend on k . This limit can easily be evaluated using asymptotic series 
to yield (for oj ^ 1) 



M 



1 ^ 



1 

3B 



1 

2B^ 



1 



3 

Ib^ 



(log, 27r + log, 2b) 



(42) 



Using series and asymptotic expansions for all the special functions present, it is routine to establish that Mj_^j_j_ « 
(48BV315)a;w'a;" for B < 1 , while for B > 1 , one finds M±^±± « uju'uj" /{3B) , the result stated in Eq. 

It must be remarked in passing that the expressions for in Eqs. (|3^), ( [A^ ) and (^2|) cannot be derived from 

the series in Eq. (C2) of pi| , principally because that series expression is divergent, and therefore erroneous. Such an 
error was introduced by an inappropriate rearrangement of individually-divergent contributing series (leading to the 
addition of infinite contributions), a mistake that is avoided by the careful technique employed here in manipulating 
the results of Appendix B in ||ll||. Notwithstanding, the numerical results for the _L^_L_L mode presented in |ll| ] 
were effectively evaluated before any series rearrangement, and therefore remain valid. 



Scaled Scattering Amplitudes for y^yy 




3 
3 



3 
3 



O 



1 
Logio[B] 

FIG. 3. The dependence of the scattering amplitudes, scaled by ujui'uo" , on the magnetic field (in units of 
Be = 4.413 X 10^^ Gauss) for the splitting modes ^^|| || (see Eq. [0) and _L-^_L_L (see Eq. [(^), for photon energies 
well below pair creation threshold. The amplitude for ||^^ || is identical to that for -L^|l || . 
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The compact analytic forms presented in Eqs. (^ij) and (^2|) represent the culmination of the w <C 1 focus here. It 
is the first time such simple forms for the scattering amplitudes involving just special functions have been calculated 
in this limit, though somewhat more convoluted, yet essentially equivalent, expressions have been put forward in 
A distinct advantage of the expressions in Eqs. (^) and (^) is the ease with which they can be accurately 
computed numerically. Their dependence on B is illustrated in Figure |^, replicatin g th e numerics of pT| and earlier 
effective Lagrangian determinations which are just as expedient since (see Eqs. [ p3| and |Q below) they involve 
just integrals of elementary functions. 

The low frequency result for 7\/tj|_,j_|| is not presented explicitly since it reproduces that for A^^^im (e.g. see 
this is due to the crossing symmetries involved. Note that while the cubic dependences of all the modes at 
low energies reflect the lack of an energy scale in this domain (i.e. such as pair threshold), the normalizations are 
dependent on the polarization mode, particularly at highly supercritical fields where the A4±^±± amplitude is highly 
suppressed. This effectively represents how the rate normalization is sensitive to the (virtual) pair creation thresholds 
for the polarization states involved in a particular splitting mode. 

To conclude this presentation focusing on the <C 1 specialization, an obvious objective is the re-derivation of 
Eqs. ( ^ll) and (^) starting with extant and well-known effective Lagrangian/proper-time (ELP) results, and thereby 
demonstrating analytically the equivalence of the S-matrix formulation in the Landau representation and Schwinger- 
type formalisms in the low energy limit. Consider first the mode J-^|| || , for which such a determination is somewhat 
involved. The starting point is the integral expression |^,^,|l^ that corresponds to the scaled scattering amplitude 
that generates the same form for the rate as in Eq. (|l2|): 

A^ELP _ ujuj'u}" f°° ds ^s/B \ ( 3 s\ coshs 3 + 2.s^ s cosh s 1 , .„^ 

-^±-1111 - \\-Ts + q) + 12 sinh^ . + 2ih^h^ / ' ^^^^ 

which has B and B ^ 1 limits matching those of Eq. (pl|). In the subsequent analysis, it is useful to manipulate 
integrations using the variable ji = 1/ B . The first step is to recognize that 1/ times the factor in curly braces in 
Eq. ( p3| ) is a perfect derivative, namely dg/ds , where g{s) = (l/4s'^) (i[scoths — l]/ds — coths/(6s) . Integration by 
parts is obviously the operative method, with the goal of retaining coth s functions explicitly, combined with powers 
of s . After some algebra, one finds that 



UlLO Ul 



1 r ds _ 



bL 



/B 

s' ^ 



2 

s coth s — 1 — g- 



(s coth s — 1) 

(44) 



results. The integral on the first line can be performed using identities 3.551.3 and 3.554.4 of |3^, yielding the Gamma 
function and polygamma functions (or equivalently generalized Riemann Zeta functions) in addition to elementary 
functions. The only subtle part pertains to the second term of this integral, namely that contributed by the — 1/ {4Bs' ) 
factor. This can be differentiated with respect to B , evaluated to yield a ip function, and then the result integrated, 
noting the behaviour as B ^ . The evaluation of the integral on the second line of Eq. (^4|) is much more involved. 
However, it has been performed before in the literature, and appears explicitly in calculations | |3^ , |38t of the one-loop 
effective Lagrangian describing refractive indices of the magnetized vacuum in QED. Hence the motivation for the 
particular partitioning of integrations chosen in Eq. (|4^). Details of the determination of this integral are found in 
Dittrich et al. [Q, and the second line of Eq. (^^ can be equated to —Stt'^ / B^ times the Lagrangian 

(see 

Eqs. (2.4) and (3.16) of [p8[ ), thereby introducing the Fi function. Collecting together the terms neatly generates 
an analytic form for A^j^m that is identical to Eq. (|4l|), so that the desired demonstration of equivalence of the 
Landau representation and effective Lagrangian forms is achieved. 

The procedure for the _L^_L_L mode is similar, though somewhat less involved. The equivalent scaled scattering 
amplitude obtained from effective Lagrangian/proper-time techniques is 



A^ELP _ Ljuj'uj" [°°ds s/B) 3 cosh 5 , 3- 4.s^ 3^^ \ 

■^^-^^--^X ~ j4iiiar7 + i;is?7"2ihd?7/- ^^^^ 

Recognizing that the factor in curly braces can be written as — (3s/4) d[cot\vs/s — l/a^j/ds -f (s^/4) d^\cot\vs — 
l/sj/ds'^, integration by parts is again indicated, with identities 3.551.3 and 3.554.4 of |33] again proving useful. 
With manipulations similar to (but simpler than: the Fi function is not involved here) those for the -L^|||| , a 
modicum of algebra leads to the derivation of Eq. ( p^ ) from Eq. (^), as desired. This equivalence is a satisfying 
indication of the verity of the Landau representation analysis in this paper. 
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IV. CONCLUSION 



This paper has provided a detailed development of the S-matrix formulation of the QED process of magnetic 
photon splitting in the Landau representation, focusing on the case of zero dispersion where photon propagation is 
collinear. The formalism in Section |l| rederives and extends the exposition of Mentzel, Berg & Wunner The 
two principal general developments offered here are an analytic reduction via the summation over the spins of the 
intermediate pair states, discussed briefly in and the analytic integration over the momenta parallel to the 
field incorporated in the electron propagators. This latter accomplishment is presented here for the first time. The 
cumulative product of these developments is a satisfyingly simple and elegant form in Eq. ( |l5|) for the scattering 
amplitude for each of the polarization modes permitted by CP invariance. These amplitudes possess products of 
generalized Laguerre polynomials that are common to QED processes in external magnetic fields, and elementary 
functions involving the photon energies and the various pair thresholds associated with the propagators. Moreover, 
the analytic forms presented consist of just triple summations over Landau level quantum numbers of the intermediate 
states, and are eminently suitable for accurate numerical computations both below and above pair creation threshold 
Lo — 2 . The applicability of these results to regimes above pair threshold is a benefit of the S-matrix expansion in the 
Landau representation that is not afforded by effective Lagrangian and proper-time calculations: while these (latter) 
Schwinger-type techniques elegantly formulate splitting rates below pair threshold, they eliminate the resonance 
structure early on in their mathematical developments, a severe limitation above oj — 2 . 

As an embellishment to these general results, specializations in two significant domains have been obtained. The 
first is for highly supercritical fields, i? » 1 , reproducing in particular the result of ||] for the -L^|||| mode, and 
deriving new results for the other two modes permitted by CP invariance in the limit of zero dispersion. The second 
group of asymptotic results are for energies w <C 1 well below pair creation threshold, where new and compact 
expressions for the scattering amplitudes have been derived in Eqs. (41) and (^) in terms of the logarithm of the 
Gamma function, its integral and their derivatives. These two domains of specialization herein have facilitated the 
first analytic demonstration of the equivalence of splitting rates obtained by the S-matrix formulation in the Landau 
representation and those derived using Schwinger-type effective Lagrangian/proper-time techniques. 
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APPENDIX A: 



Here the remainders that appear in the series representation in Eq. ( p8| ) for the a; <C 1 specializations to the 
splitting amplitudes for polarization modes -L^|||| and _L^_L_L are presented: 
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and 
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7^i^ii(fc, B) = - 
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APPENDIX B: 



In this Appendix, various useful properties of the ri(a;) function, the integral of the logarithm of the Gamma 
function, that are needed in the w ^ 1 specializations are stated. Given the definition of Fi in Eq. (^), it is 
elementary to establish, using 44.1.2 of ]B§], that 



n=0 



(x + n)log^{x + n) = loggri(l 



fc)-log,ri(x) 



(Bl) 



Taking successive derivatives with respect to x , one quickly arrives at well-known finite series representations of T(x) 
and its logarithmic derivative ip{x) ; see |33| , ^ for discussions of these functions and their series representations. 

An asymptotic series representation for the Fi function for large arguments is useful, and can be derived with the 
aid of the following series representation (see result 8.343.2 of [53) for the logarithm of the Gamma function: 



1 1 1 ^ 

loge r(2;) = ~ 2 ) ^^^'^ a; - X + 2 loge 271- + 2 XI 



(m + 1)(to + 2) 



({m. + 1, x) 
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(B2) 



from which Stirling's asymptotic expansion can be derived. Here, CC*^: is tti^ generalized Riemann Zeta function, 
defined in 9.511 and 9.521.1 of [Q. The integration of this series is effected using the identity C'("^j t) = —({m+1, t) , 
and is mostly uneventful. However, the treatment of the m = 1 term in the summation is somewhat more subtle, due 
to the singular nature of C(Oj^) j ^-iid requires taking the limit m — > l"*" , assuming m to be a continuous variable. 
Then result 8.362.1 of pSl comes in handy, and the series identity 
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(B3) 



follows, where 7e = — V'(l) ~ 0.5772 is Euler's constant. This series, which adequately substitutes for an asymptotic 
representation, can be used very effectively for numerical evaluations for all a; > 1 . For the range < x < 1 , this series 
also effects accurate evaluation of logg Fi {x) via use of the recurrence relation logg Fi (x) — logg Fi (1 + x) — x logg x , 
an identity derivable from Eq. (40) with the aid of 6.441.3 in [Q. For large x , it then follows that 
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where 



Li 



ill V C(m, 1)-1 
12 ' 6 2 (m + 1)(to + 2) 



fc(fc + 1) 
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(B5) 



with numerical value Li « 0.24875 . This is just the constant appearing in the magnetized vacuum polarization 
analyses of [ ^7||3^ , where the Raabe integral form for it can be found. The second definition of Li in Eq. (B5) can 
be obtained by setting x = in Eq. (Bl), and is a result noted by El|. 
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